Models of magnetized neutron star 
atmospheres 



V. Suleimanov 

Institut fur Astronomic und Astrophysik, Kepler Center for Astro and Particle 
Physics, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany; 
Kazan State University, Kremlevskaja str., 18, Kazan 420008, Russia 

A.Y. Potekhin 

loffe Physical-Technical Institute, Polytekhnicheskaya str., 26, St. Petersburg 

194021, Russia 

K. Werner 

Institut fur Astronomie und Astrophysik, Kepler Center for Astro and Particle 
Physics, Universitat Tubingen, Sand 1, 72076 Tubingen, Germany 



Abstract 

We present a new computer code for modeling magnetized neutron star atmo- 
spheres in a wide range of magnetic fields (10 12 — 10 15 G) and effective tempera- 
tures (3 x 10 5 — 10 7 K). The atmosphere is assumed to consist either of fully ionized 
electron-ion plasmas or of partially ionized hydrogen. Vacuum resonance and partial 
mode conversion are taken into account. Any inclination of the magnetic field rela- 
tive to the stellar surface is allowed. We use modern opacities of fully or partially 
ionized plasmas in strong magnetic fields and solve the coupled radiative transfer 
equations for the normal electromagnetic modes in the plasma. Using this code, we 
study the possibilities to explain the soft X-ray spectra of isolated neutron stars 
by different atmosphere models. In particular, the outgoing spectrum using the 
"sandwich" model (thin atmosphere with a hydrogen layer above a helium layer) is 
constructed. Thin partially ionized hydrogen atmospheres with vacuum polarization 
are shown to be able to improve our understanding of the observed spectrum of the 
nearby isolated neutron star RBS 1223 (RX J1308.8+2127). 
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1 Introduction 



At present several classes of neutron stars (NSs) with strong magnetic fie 



are known. They include X-ray dim isolated NSs (XD INSs. lHaberl 



2007) 



20041 ) 



sec 



central compact objects (CCOs) in supernova remnants (jPavlov et al. 
anomalous X-ray pulsars and soft-gamma repeaters (AXPs an d SGRs, 
reviews bv iKaspil 120071 ; iMereghetti et alJlioOTl lMereghettill2008h . The NSs in 
these classes have superstrong magnetic fields up to B > 10 14 G (SGR and 
AXP) and B ~ a few x 10 13 G in XDINSs, evaluated from period changes and 
from absorption features in the observed spectra, if they are interpreted as ion 
cyclotron lines (see reviews by lHaberll 120071 ; Ivan Kerkwijk fc Kaplan! 120071 ). 



These NSs are sufficiently hot (T e g ~ 10 6 — 10 7 K) to be observed as ther- 
mal soft X-ray sources. Some of the XDINSs and CCOs have one or more 
absorption fea tures in their X-ray spectrum at the energies 0.2 - 0.8 keV 



(lHaberll . 120071 ) , an d the central energi e s of these features appear to be har 



20071 : Ivan Kerkwiik fc Kaplanl . 120071 : lHaberll . l20oj )7The optical/ultraviolet 



appj 

moni c ally spaced dSanwal et all. 120021; iBignami et all 120031 ; ISchwope et all 



fluxes of the known XDIN optical counterparts a re a few times large r than 



Kaplan et al.. 


2003 


:lMotch et al.. 


2003; 



Mignani et all 12007( 1 . 



2001, 2003 



The XD INs are nearby object s, and parallaxes of some of them have been 
measured (IKaplan et all l2002al ). Therefore, they give a good possibility to 
measure the NS radii, yiel ding useful information on the equation of stat e 
(EOS) for the NS inner core (ITrumper et all 12004 ; llattimer &: Prakashl . l2007l ). 
For a sufficiently accurate evaluation of NS radii, a good model of the NS sur- 
face radiation is necessary for the observed X-ray spectra fitting. 



Structures and emergent spectra of NS atmospheres with str ong (B > 10 



12 



G) magnetic fields haye been modeled by many scientific groups (IShibanov et al 



1992 



Raiagopal etafl Il997l : lOzell . l200ll ; IHo fc lail . l200ll . 120031 . 12004 Ivan Adelsberg fc lai 



2006). Methods of fully ionized model atmospheres modeling (see, e.g.. IZavlin 
2009 for references) and partially ionized hydrogen atmospheres modeling 
flPotekhin etHL I2004J : IHo fc Lai 12004 IHo et all l2008f ) were developed. Mid- 
Z element a t mosp heres for strongly magnetized NSs have also been modeled 
(jMori fc Hoi . I2007I ). The effe ct of the vacuum polarization on magnetized N S 
atmospheres was studied by iPavlov fc Gnedinl (119841 ) ; llai fc Hoi (120021 . 120031 ) , 
and model atmospheres wi th partial mode conversion due to vacuum p olariza- 
tion have been computed (IHo fc lail . 120031 ; Ivan Adelsberg fc lail . 120061 ) . Some 
of the XDINSs might have a "thin" atmosphere above the condensed surface, 
which could be optically thick to low-energy photons and optically thin to 
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high-energy photons (IMotch et al 



2003 



Ho et al 



20071). 



Here we present a new computer code, which allows to compute magnetized 
NS model atmospheres consisting of partially ionized hydrogen, taking into 
account the vacuum polarization effect together with partial mode conversion 
and an arbitrary magnetic field inclination. Compared to previous codes, other 
methods were used for the temperature correction procedure and the solution 
of the radiation transfer equation. We also present some new results, which 
were obtained by using this code. 



2 Method of atmosphere structure calculations 



We computed model atmospheres of hot, magnetized NSs subject to the 
constraints of hydrostatic and radiative equilibrium assuming planar geometry. 
There are two versions of the code: the first one is designed for the magnetic 
field B perpendicular to the surface, and in the second version, the angle 6 b 
between B and n is arbitrary, and calculations are more expensive. 

The model atmosphere structure for an NS with effective temperature T e g, 
surface gravity g, magnetic field B, and given chemica l composition is de- 



scribed by a set of equations (e.g., IShibanov et al.lll992l ; iHo fc Laill200ll and 



references therein), which are written below for the simplest first case , where 



9b = 0. In the second case (for arbitrary 9b, e.g., Zavlin et al. 19951 ). there 



is an additional dependence of the absorption {k l v ) and scattering (al) coeffi- 
cients on the azimuthal angle accordingly, there are additional integrations 
over ip in Eqs. (3), (6), (8), and (9) in the second case. 

The hydrostatic equilibrium equation reads 
dP e 



dm 



9 ~ fi'rad, 



where 



GM NS 

9 = 2 r— (2) 

-^nsv 1 — Rs/Rns 

2 oo +1 



fed = — £ / dv f (** + oi) ii PM dfx, (3) 
c i=1 J J 



and Rs = 2GMns/c 2 is the Schwarzschild radius. is the specific intensity 
in mode i, fi — cos 9, where 9 is the angle between the surface normal and 
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the radiation propagation direction, P g is the gas pressure, and the column 
density m is determined as 

dm = —p dz . (4) 



The variable p denotes the gas density and z is the vertical distance. 
The radiation transfer equations for the two modes are 

where the source function is 

and B v is the blackbody (Planck) intensity. The optical depth r* is defined as 
drt = (K + <) dm. (7) 

Of course, opacities k v and a l v depend on p. These equations have to be com- 
pleted by the energy balance equation 

oo 2 +1 

/ duJ2 j ((K + <)I l M - VM) dp = (8) 

o i=1 -i 

with emissivity 

2 + 1 



VM = K\ + £ / <tf'fo A*0 dp\ (9) 

j'=i_i 



Equations (1) - (9) must be completed by the EOS and the charge and 
particle conservation laws. In the code two different cases of these laws are 
considered. 

In the first (simplest) fully ionized atmosphere is calculated. There- 

fore, the EOS is the ideal gas law 

P g = n tot kT, (10) 



4 



where n tot is the number density of all particles. Opac i ties a re calculated in 
the same way as in the paper by Ivan Adelsberg fc Lail (l2006h (see references 
therein for the background theory and more sophisticated approaches). The 
vacuum polarization effect is taken into account following the same work. 



The second considered case is a partially ionized hydrogen atmosphere. In 
this case t he EOS and the corresponding opa cities are taken from tables cal- 
culated by iPotekhin fc Chabrierl (120031 . 120041) . The normal mod e polarization 
vectors are taken from the calculations by lPotekhin et al.l (12004 ). The vacuum 
polarization effect is also included. 



We use a logarithmically equidistant photon energy set in our computations 
in the range 0.001 - 20 keV with additional points near each ion cyclotron 
resonance. If the vacuum resonance is taken into consideration, then another 
photon energy gri d is used, which is constructed using the "equal grid" method 
(IHo fc Lail . 120031 ). This energy grid is recalculated after every iteration. 



The radiation transfer equation (5) is solved on a set of 40 polar angles 9 
(and additionally 6 azimuthal angle s <p in the case of inclin ed magnetic field) 
by the short characteristic method (jOlson fc Kunaszl . 119871 ). We use the con- 
ventional condition (no external radiation) at the outer boundary. At the inner 
boundary we take the incoming specific intensities equal to blackbody radia- 
tion corresponding to the temperature of the last atmosphere point. In the case 
of the thin atmosphere the temperature at the inner boundary is considered 
as the temperature of a condensed NS surface. The co de allows one to take 



i nto a ccount the partial mode conversion according to Ivan Adelsberg fc Lai 
fl2006h . 



The solution of the radiative transfer equation (5) is checked for the energy 
balance equation (8). Temperature corrections are then evaluated using three 
different procedures (IKuruczl . Il970l ). modified to deal with strong magnetic 
fields. A nonmagnetic versi on of this code was describ ed, for example, by 
Ibragimov et al.l (120031 ) and ISuleimanov fc Wernerl (120071 ) . 



The iteration procedure is repeated until the relative flux error is smaller 
than 1% and the relative flux derivative error is smaller than 0.01%. 



Our method of calculation ha s been tested b y a comparison to models 



for magnetized NS atmospheres (IShibanov et all Il992l ; iPavlov et al.l . 11994 



Ho fc Lai l200ll ; lOzell . l200ll : IHo fc Lai l2003h . Model atmos pheres wi t h par 



tially ionized hydrogen are compared to models computed by lHo et all (120071 ). 
We have found that our models are in a good agreement with these calcula- 
tions. 
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Fig. 1. Emergent spectra of the thin fully ionized hydrogen atmospheres above a 
solid surface with T e g = 10 6 K and B = 10 14 G in comparison to the semi-infinite 
atmosphere (dashed curve). The spectra of atmospheres with surface densities S=l 
(solid curve) and 100 g cm -2 (dash-dotted curve) together with the corresponding 
blackbody (dotted curve) are shown. 

3 Results 



Here we present some new results, obtained by using of the developed code. 
In all calculations below we use the same surface gravity, log g = 14.3. 

One of the problems relate d to AXPs is the lack of any absorption feature at 
the proton cyclotron energy. iHo fc Lail (120031 ) suggested that a possible solu- 
tion of this problem is the suppression of the cyclotron absorption feature due 
to the vacuum polarization. We demonstrate that this absorption line is also 
further reduced in a thin atmosphere. Figure 1 demonstrates emergent spectra 
of thin atmospheres without allowance for the vacuum polarization effect. The 
absorption feature disappears with decreasing the atmosphere surface density 
E. Equivalent widths of the absorption feature are ~ 670 eV for the semi- 
infinite atmosphere, m 190 eV for the X = 100 g cm -2 slab and ~ 5 eV for 
the E = 1 g cm -2 slab. This decrease of the equivalent width corresponds to 
the decrease of the optical depth according to the curve of growth. 

Various hypotheses were considered for an explanation of the two harmoni- 



cally spaced absorption features in CCOs (see lMori fc Holl2007l and references 
therein). Here we suggest another one, which we name "sandwich atmosphere" . 
A thin, chemically layered atmosphere above a condensed NS surface can arise 
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Fig. 2. Emergent spectrum (top panel) and temperature structure (bottom panel) 
of the "sandwich" model atmosphere above a solid surface with T e g = 10 6 K and 
B = 10 14 G (solid curve) in comparison with the thin fully ionized hydrogen at- 
mosphere (dashed curve) with the same parameters. The surface densities of both 
model atmospheres are 100 g cm -2 , in the "sandwich" model the H slab has 25 g 
cm -2 surface density and the He slab has 75 g cm -2 . The corresponding blackbody 
spectrum (dotted curve) is also shown in the top panel. 
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from accretion of interstellar gas with cosmic chemical composition. In this 
case, hydrogen and helium quickly separate due to the high gravity. In this 
"sandwich atmosphere" a layer of hydrogen is located above a helium slab, and 
the emergent spectrum has two absorption features, corresponding to proton 
and a-particle cyclotron energies. In Fig. 2 the emergent spectrum for one of 
such models together with the corresponding temperature structure are shown. 
Fully ionized hydrogen and helium layers are considered in this model. The 
emission feature at the helium absorption line arises due to a local temper- 
ature bump at the boundary between the helium and hydrogen layers. This 
temperature increase compensates a lower number density in the hydrogen 
layer to avoid a gas pressure jump. Of course, this temperature jump should 
be smoothed by hydrogen-helium mixing and thermal conduction flux at the 
boundary between layers. Here we used the sharp boundary and did not take 
into account thermal conductivity. Results can also change quantitatively, if 
partially ionized hydrogen and helium layers are considered. 



Most of the X DINSs have m agnetic fields B > 10 13 G and color temper- 
atures ~ 10 6 K ( IHaberll . 120071 ). Hydrogen model atmospheres are partially 
ionized under these conditions and the vacuum polarization effect is also sig- 
nificant. Here we present first results of modeling of partially ionized hydrogen 
atmospheres with the vacuum polarization effect and partial mode conversion 
using our radiative transfer code. Previously part ial mode conversion was con - 
sidered only for fully ionized atmosphere models (Ivan Adelsberg fc Lail . 120061 ) . 
In Fig. 3 we compare spectra and temperature structures of the partially ion- 
ized hydrogen model atmospheres with and without the partial mode conver- 
sion effect. When the X-mode (having smaller opacity) partially converts to 
the O-mode in the surface layers of the atmosphere, the energy absorbed by 
the O-mode heats these upper layers. As a result, the emergent spectra are 
closer to the blackbody. The absorption feature to the right of the proton cy- 
clotron line arises mainly due to the transition from the ground state to the 
first excited state of the H atom; it has a maximum at 350 eV and is strongly 
broadened at lower energie s because of the atomic motion effect (cf. Fig. 5 in 



Potekhin fc Chabrierll2004j and related discussion there). 



The effective temperature and magnetic field strengths are not uniform 
over the NS surface, and generally the magnetic field is not perpendicular to 
the surface. Therefore, for comparison to observations, it will be necessary 



to integrate the loc al model spectra over the NS surface (see IZavlin et al 



19951 ; iHo et al.ll2008l ). and to compute model atmospheres with inclined mag- 



netic field. This possibility is included in our code. For example, Fig. 4 shows 
spectra and temperature structures of partially ionized hydrogen model at- 
mospheres with magnetic field perpendicular and parallel to the NS surface. 
The vacuum polarization effect is taken into account without mode conver- 
sion. The difference between emergent spectra is not si gnificant (see also t he 



first models with different magnetic field inclinations in lKaminker et al.lll982 
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Fig. 3. Emergent spectra and temperature structures of the partially ionized hy- 
drogen model atmospheres with T e g = 1.5 x 10 6 K with (solid curves) and without 
(dashed curves) vacuum polarization effect (with the partial mode conversion) . The 
magnetic field strength is B = 4 x 10 13 G. The corresponding blackbody spectrum 
is also shown in the upper panel (dotted curve). 
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Fig. 4. Emergent spectra and temperature structures of partially ionized hydro- 
gen model atmospheres with T e g = 1.5 x 10 6 K with inclinations of the magnetic 
field (B = 4 x 10 13 G) to the surface normal 6b equal 0° (solid curves) and 90° 
(dashed curves). The vacuum polarization effect is not included. The corresponding 
blackbody spectrum is also shown in the upper panel (dotted curve). 
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Fig. 5. Emergent spectra and temperature structures of the partially ionized hydro- 
gen model atmospheres with T e g- = 1.2 x 10 6 K with vacuum polarization effect and 
partial mode conversion shown for various surface densities £ (solid curves - 1 g 
cm -2 , dashed curves - 10 g cm -2 , dash-dotted curves - semi-infinite atmosphere). 
The magnetic field strength is B = 4 x 10 13 G. The corresponding blackbody spec- 
trum is also shown in the upper panel (dotted curve). 
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Shibanov et al.l 1 19921) , but it can be more significant at other atmospheric 
parameters (IHo et al.l . 120081 ) . 



Ho et al.l (120071 ) demonstrated that a single partially ionized thin hydrogen 
atmosphere can explain the problem of the observed optical/UV flux excess 
over the blackbody extrapolation from X-ray to the optical range in the case 
of brightest isolated NS RX J1856.4-3754: the model fits well both the ob- 
served optical flux and the X-ray spe ctrum. RX J1856.4— 3754 ha s very low 
pulsed fraction of radiation (« 1.2%, iTiengo &: Mereghettil 120071 ) . therefore 
it is possible to fit the radiation of this star by a single model atmosphere 
(although the small pulsed fraction can also be explained by a small value of 
an angle between the magnetic and rotation axes or between the rotation axis 
and the lin e of sight). Other X DINSs have larger pulse fractions, up to 18% 
(RBS 1223, lHaberl et al.l 120041 ) . In this case the temperature is certainly not 
uniform across the NS surface, and the excess optical flux ca n be explained 



by the radiation from cool surface parts (jSchwope et al.l . 120051 ) 



We now investigate properties of partially ionized hydrogen models, which 
can be applied to the RBS 1223 atmosphere. The color temperature of this 
star, found from X-ray spectra fi t ting, is close to 10 6 K, with magnetic field 
B ~ 4 x 10 13 G ( jSchwope et al.l . 120071 ) . In particular, we study the optical 
flux excess in comparison to the X-ray fitted blackbody flux in this kind of 
models in order to explore which part of the observed optical excess can be 
explained by the atmosphere effect. This optical excess is illustrated in the 
top panel of Fig. 6. For this aim we have calculated two sets of models with 
vacuum polarization and partial mode conversion. The models in the first 
set have effective temperatures T c s = 10 6 K and the models of second one 
have effective temperatures T e g = 1.2 x 10 6 K. In both sets B = 4 x 10 13 
G, and models with surface densities S = 1, 3, 10, 30, 100 and 10 5 (semi- 
infinite model) g cm -2 are computed. In Fig. 5 we show emergent spectra 
and temperature structures for some models from the second set. Clearly, the 
X-ray spectra of the models with E < 10 g cm -2 are close to a blackbody and, 
therefore, better fit the observed X-ray spectrum. Let us re mark that there 



is a s ignificant absorption feature in the RBS 1223 spectrum (jSchwope et al. 



20071 ). which cannot be fitted by the simple blackbody spectrum. Here we do 
not try to fit the observed spectrum of RBS 1223 by the single temperature 
model: clearly, such fitting would require the integration of the local spectra 
over the neutron star surface ass uming some non-uniform temperat ure and a 



magnetic field distributions (see Zavlin et al. 1995 ; Ho et al.l 20081 ). Instead 



we are trying to understand, what kind of atmosphere (thin or semi-infinite) 
be more appropriate for this modeling. 



In Fig. 6 (bottom panel) we show the ratio of the model atmosphere flux to 
the X-ray fitted (in the 0.4-1 keV range) blackbody flux in the optical band 
at A = 5150 Adepending on E for both sets. The observed ratio is about 5 
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Fig. 6. Top panel: Emergent spectrum of the partially ionized hydrogen model at- 
mosphere of neutron stars with T c g = 10 6 K with vacuum polarization effect, par- 
tial mode conversion and magnetic field strength B = 4 x 10 13 G. The blackbody 
spectrum fitted to the maximum of the spectral distribution is also shown (dashed 
curve) . At the optical band the model atmosphere flux is a few times larger than the 
blackbody flux. Bottom panel: Ratios of the model atmosphere flux to the black- 
body (X-ray fitted) flux at the optical band depending on the model atmosphere 
thickness (surface density S) for the models with T e g- = 10 6 K and T e fj = 1.2 x 10 6 
K. 
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( jKaplan et al.l . l2002bl ). in agreement with the semi-infinite atmosphere mod- 
els. However, the observed blackbody like X-ray spectrum agrees with the thin 
atmosphere models, for which this ratio is close to 1. Therefore, the observed 
optical excess cannot be explained by the thin atmosphere model alone; in- 
stead, it can arise due to a nonuniform surface tem perature distribution , in 
agreement with the RBS 1223 light curve modeling (ISchwope et al.l . 120051 ) . 



4 Conclusions 



In this work we present a new code, which can model fully ionized and par- 
tially ionized hydrogen atmospheres in a wide range of effective temperatures 
(3 x 10 5 - 10 7 K) and magnetic fields (10 12 - 10 15 G), with any inclination 
of the magnetic field to the stellar surface. The vacuum polarization effect 
with partial mode conversion is taken into consideration. Calculated emer- 
gent spectra and temperature structures of the model atmospheres agree with 
previously published ones. 

We have studied the properties of thin atmospheres above condensed NS 
surfaces. We demonstrated that the proton cyclotron absorption line disap- 
pears in the thin hydrogen model atmospheres. A new thin "sandwich" model 
atmosphere (hydrogen layer above helium layer) is proposed to explain the 
occurrence of two absorption features in the observed X-ray spectra of some 
isolated NSs. 

A set of of partially ionized hydrogen model atmospheres with vacuum 
polarization and partial mode conversion with parameters (effective tempera- 
ture and the magnetic field strength) close to the probable parameters of the 
isolated NS RBS 1223 were calculated. We analysed the optical excess (rela- 
tive to the X-ray fitted blackbody flux) in the model spectra and found that 
the optical flux excess ~ 5 for the semi-infinite model atmospheres decreases 
down to 1 with decreasing surface density £ of the atmosphere. Spectra of 
thin model atmospheres are closer to the observed RBS 1223 X-ray spectrum, 
therefore we conclude that the observed optical excess should be explained by 
nonuniform surface temperature distribution. 
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